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ABSTRACT 

The encroachment of woody plants in grasslands across the Western U.S. will affect soil 
water availability by altering the contributions of evaporation (E) and transpiration (T) to total 
evapotranspiration (ET). To study this phenomenon, a network of flux stations is in place to 
measure ET in grass- and shrub-dominated ecosystems throughout the Western U.S. A method 
is described and tested here to partition the daily measurements of ET into E and T based on 
diurnal surface temperature variations of the soil and standard energy balance theory. The 
difference between the mid-afternoon and pre-dawn soil surface temperature, termed Apparent 
Thermal Inertia (Ia), was used to identify days when E was negligible, and thus, ET=T. For 
other days, a three-step procedure based on energy balance equations was used to estimate the 
contributions of daily E and T to total daily ET. The method was tested at Walnut Gulch 
Experimental Watershed in southeast Arizona based on Bowen ratio estimates of ET and 
continuous measurements of surface temperature with an infrared thermometer (IRT) from 2004- 
2005, and a second dataset of Bowen ratio, IRT and stem- flow gage measurements in 2003. 
Results showed that reasonable estimates of daily T were obtained for a multi-year period with 
ease of operation and minimal cost. With known season-long daily T, E and ET, it is possible to 
determine the soil water availability associated with grass- and shrub-dominated sites and better 
understand the hydrologic impact of regional woody plant encroachment. 
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POPULAR SUMMARY 

In the semiarid and arid regions of the western United States, woody plants in the form of shrubs 
have begun to encroach upon traditional grasslands. It is known that shrub-dominated vegetation 
has a different water demand than herbaceous vegetation, which affects the balance between soil 
evaporation (E) and plant transpiration (T). This is turn impacts soil water availability, which is 
the driving force behind this region’s biologic, hydrologic, and socioeconomic processes. 
However, accurate partitioning of water loss between E and T is one of the most important 
ecohydrological challenges in understanding vegetation dynamics in dryland environments. 
With better seasonal information about E and T, it should be possible to better interpret the 
biological and hydrological impacts of woody shrub encroachment and plan for the subsequent 
socioeconomic consequences related to this land use change. 

This paper describes a method to partition the daily measurements of evapotranspiration or ET 
(total evaporation plus transpiration) which are available from a network of flux stations 
throughout the western U.S. into E and T based on diurnal surface temperature variations of the 
soil and standard energy balance theory. Results showed that reasonable estimates of daily T 
were obtained for a multi-year period with ease of operation and minimal cost. With known 
values for season-long daily T, E, and ET such as obtained through this study, it is possible to 
determine the soil water availability associated with grass- and shrub-dominated sites and better 
understand the hydrologic impact of regional woody plant encroachment. 
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SIGNIFICANT FINDINGS STATEMENT 

Background: In the semiarid and arid regions of the western United States, woody plants in the 
form of shrubs have begun to encroach upon traditional grasslands. It is known that shrub- 
dominated vegetation has a different water demand than herbaceous vegetation, which affects the 
balance between soil evaporation (E) and plant transpiration (T). This is turn impacts soil water 
availability, which is the driving force behind this region’s biologic, hydrologic, and 
socioeconomic processes. However, accurate partitioning of water loss between E and T is one 
of the most important ecohydrological challenges in understanding vegetation dynamics in 
dryland environments. With better seasonal information about E and T, it should be possible to 
better interpret the biological and hydrological impacts of woody shrub encroachment and plan 
for the subsequent socioeconomic consequences related to this land use change. 

Approach: To address this issue, a method was developed to partition the daily measurements 
of evapotranspiration or ET (total evaporation plus transpiration) which are available from a 
network of flux stations throughout the western U.S. into E and T based on diurnal surface 
temperature variations of the soil and standard energy balance theory. The approach involved 
the addition of infrared thermometers to standard flux stations, and using a thermal inertia (I A ) 
calculation (difference between midday and predawn surface temperatures) to identify days 
when evaporation from the soil was negligible. Results showed that reasonable estimates of 
daily T were obtained for a multi-year period with ease of operation and minimal cost. 

Significance: The main contribution of this work was a new method to partition on-site 

measurements of daily total evapotranspiration (ET D ) into daily E and T based on the low-cost 
addition of an infrared thermometer to existing eddy covariance and Bowen ratio stations. 
Reasonable estimates of To were obtained for clear-sky days when Ed was found to be zero 
based on the magnitude of the apparent thermal inertia (I A ). The approach is based on 
instrumentation that can be maintained in place continuously for years with no more expertise 
and effort than is already required for deployment of energy flux stations. 
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ABSTRACT 

The encroachment of woody plants in grasslands across the Western U.S. will affect 
soil water availability by altering the contributions of evaporation (E) and transpiration 
(T) to total evapotranspiration (ET). To study this phenomenon, a network of flux 
stations is in place to measure ET in grass- and shrub-dominated ecosystems throughout 
the Western U.S. A method is described and tested here to partition the daily 
measurements of ET into E and T based on diurnal surface temperature variations of the 
soil and standard energy balance theory. The difference between the mid-afternoon and 
pre-dawn soil surface temperature, termed Apparent Thermal Inertia (Ia), was used to 
identify days when E was negligible, and thus, ET=T. For other days, a three-step 
procedure based on energy balance equations was used to estimate the contributions of 
daily E and T to total daily ET. The method was tested at Walnut Gulch Experimental 
Watershed in southeast Arizona based on Bowen ratio estimates of ET and continuous 
measurements of surface temperature with an infrared thermometer (IRT) from 2004- 
2005, and a second dataset of Bowen ratio, IRT and stem-flow gage measurements in 
2003. Results showed that reasonable estimates of daily T were obtained for a multi-year 
period with ease of operation and minimal cost. With known season-long daily T, E and 
ET, it is possible to determine the soil water availability associated with grass- and shrub- 
dominated sites and better understand the hydrologic impact of regional woody plant 
encroachment. 


INTRODUCTION 

Encroachment of woody plants in grasslands has become a common phenomenon 
across the Western U.S. over the past 150 years (Schlesinger et al., 1990; Van Auken, 
2000). This is of particular ecologic interest because these grassland communities have 
not been invaded by nonnative species, but rather indigenous species have increased 
because of changes in local abiotic or biotic conditions. Furthermore, there is serious 
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concern that this trend will continue and the change in vegetation composition and 
structure may be irreversible (Peters et al., 2004). 

The impact of this transition on ecosystem function is still unclear, resulting in a 
multitude of studies contrasting the biological processes in grass- and shrub-dominated 
sites. These studies have addressed biodiversity (Murphy and Weiss, 1992), vegetation 
distribution (Guswa et al. 2002; Smith et al. 1995), plant productivity (Lauenroth and 
Sala, 1 992), and plant water use efficiency (Emmerich and Verdugo, this issue (a)). The 
physical sciences have received similar consideration with studies comparing grassland 
and shrubland runoff and erosion patterns (Abrahams et al., 1995; Bhark and Small, 

2003; Wilcox et al., 2003), streamflow (Wilcox, 2002), soil fertility (Schlesinger et al., 
1990; Kieft et al., 1998), energy balance (Small and Kurc, 2003), and soil moisture and 
texture (Scott et al., 2000; Femandez-Illescas et al., 2001; Yao et al., in press). Results 
have shown that the grass-to-shrub transition could have far-reaching impacts on all 
aspects of ecosystem function ranging from soil microfauna (Schlesinger et al., 1990) to 
regional rainfall patterns (Taylor, 2000). This, in turn, could have socioeconomic 
impacts related to land use change as well as climate change (Houghton et al., 1999). 

In semiarid and arid regions, these studies have focused on soil water availability, 
which is the driving force behind ecosystem biologic, hydrologic and socioeconomic 
processes. Shrub-dominated vegetation has a different water demand from that of 
herbaceous vegetation, manifesting in the water loss from evapotranspiration (ET). 

There is evidence that woody plant encroachment may not impact the total ET (Kurc and 
Small, 2004; Dugas et al., 1996; Phillips, 1992), but it can alter the relative contributions 
of soil evaporation (E) and plant transpiration (T) to ET. In turn, these shifts in E versus 
T related to vegetation change can impact net ecosystem production and carbon cycling. 
This has important implications for resource management strategies and other surface 
manipulations in dryland ecosystems associated with intensification of land use and 
climate change (Loik et al., 2004). In a landmark analysis of vegetation dynamics in 
drylands, Huxman et al. (2005) identified the partitioning of E and T as one of the most 
important ecohydrological challenges in understanding vegetation dynamics in drylands. 
That is, with seasonal information about E and T under current conditions, it may be 
possible to better interpret the biological and hydrological impacts of shrub encroachment 
and plan for the socioeconomic consequences. 

Reynolds et al. (2000) summarized the results of multiple studies in aridlands and 
reported the very different measurements of the percentage of total ET attributed to T, 
varying from as little as 7% to as much as 80% in arid and semiarid ecosystems in the 
American Southwest (Arizona, Nevada and New Mexico). Nonetheless, these studies 
made the first attempts to interpret the T/ET ratio in terms of vegetation spatial patterns 
and energy balance (Tuzet et al., 1997), precipitation (Taylor, 2000; Loik et al., 2004), 
carbon dioxide exchange (Scott et al., in press), geomorphology (Smith et al., 1995) and 
plant community (Schlesinger et al., 1987). 

These studies were carried out over relatively short time periods and did not 
represent the variability of ET associated with seasonal variability in rainfall. Lane et al. 
(1983) and Reynolds et al. (2000) used model simulations to better understand the T/ET 
ratio of desert sites over longer time periods (nine- and 1 00-year periods, respectively), 
and still reported values of T/ET varying from 1-58% for grassland and 6-60% (or 15- 
37%) for shrub-dominated sites. However, with these longer data sets, they could 
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interpret the effects of seasonal precipitation on plant transpiration and net production. 
Reynolds et al. (2000), in particular, emphasized the link between T/ET and rainfall 
patterns that explained the great year-to-year variability in results for grass- and shrub- 
dominated sites. Furthermore, they were able to simulate environmental conditions 
before and after a “treatment” and to predict the hydrologic consequences of various land 
uses and management practices. Such simulation models and others (e.g. Laio et al., 
2001) offer a handy tool for partitioning ET and evaluating trends and relative variations, 
however the accuracy of absolute results is uncertain. Guswa et al. (2002) reported that 
ET estimated with a simple bucket-filling model and a more complex model at one 
location differed by 50% and the T/ET ratio varied from 55% to 67% over a season. 

The goal of this work is to offer a new way to measure, rather than simulate, the 
T/ET over long time periods with ease of operation and minimal cost. An approach is 
proposed in which daily ET, measured with a conventional eddy covariance or Bowen 
ratio technique, is partitioned into E and T using coincident measurements of diurnal soil 
surface temperature and basic energy balance theory. Instrumentation for these 
measurements can be maintained in place continuously for years, as demonstrated in this 
and other studies. Then, with known season-long daily T, E and ET, it is possible to 
determine the soil water availability associated with grass- and shrub-dominated sites and 
better understand the hydrologic impact of regional woody plant encroachment. 

APPROACH 

This approach is based on the assumption that conventional eddy covariance or 
Bowen ratio instrumentation is in place at a site, with coincident measurements of soil 
surface temperature, making measurements throughout the day at the commonly used 20- 
or 30-minute time interval. In this approach, the difference between the mid-afternoon 
and pre-dawn soil surface temperature, termed Apparent Thermal Inertia (U), was used to 
identify days when E was negligible. It is demonstrated herein that when Ia reached a 
seasonal maximum, E approached zero. With this set of measurements at a given site, 
daily ET (ETd) can be measured; dates for which daily E approaches zero (Ed~ 0) can be 
identified; and daily T (T D ) can be equal to ETd, where To = ETd- 

It is shown here that this approach can provide reliable estimates of To for 
cloudfree days when Ed~0. For days not meeting these two criteria, a three-step 
procedure based on energy balance equations was used to estimate the contributions of 
Ed and Td to total ETd- First, for clear-sky dates when ET D ~0, a value of aerodynamic 
resistance was computed based on the assumption that all available energy was converted 
to sensible heat (see details below). Second, the thus-calibrated energy balance equations 
were used to compute an E Index (El) defined as the ratio of E to potential E (E p ) at 
midday, and actual E was computed as E=E P *EI. Third, Td was computed for all “other” 
days (when Ed 5$ and/or skies were not cloudfree) by Td= ETd- Ed- 

With this set of measurements and equations, it was possible to partition 
measurements of ETd into To and Ed for all days in the growing season. Background 
information about thermal inertia and energy balance equations is given in the following 
subsections. 

Apparent Thermal Inertia 
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By definition, soil thermal inertia (I) represents the ability of soil to conduct and 

store heat, where 

I = (kpc) 1/2 [J nf 2 K' 1 s' 1/2 ]. (1) 

In Eq. (1), k = thermal conductivity [W m' K' ] ; p = density [kg m' ]; and c = heat 
capacity [J kg' 1 K' 1 ] . Like I, apparent thermal inertia (I A ) also represents the resistance 
of soil to temperature change. However, it is derived instead from the difference between 
mid-afternoon and pre-dawn surface (or soil) temperatures, where 

Ia — (ts2pm “ ts5am) [ C] . (2) 

The terms t S 2 P m and t S 5 am represent soil surface temperatures measured with a down- 
looking infrared thermometer (IRT) at times 2:00 pm and 5:00 am, respectively. 

In early studies, I A was loosely related to regional soil moisture (Kahle et al., 

1987; Pratt and Ellyett, 1979). Though introduced in the early 1980s based on satellite 
images of surface temperature (Price, 1977), it was not easily interpreted over a 
heterogeneous terrain (Price, 1985). That is, I A responds to changes in soil moisture and 
mineralogy, but it is also highly sensitive to changes in incoming solar radiation, as well 
as wind speed, air temperature and vapor pressure. In this application, these fundamental 
limitations in application of I A are overcome by 1) computing I A at one site and 
interpreting the signal over time rather than space, and 2) combining I A with an on-site 
measurements of the surface (in this case, ET) to account for atmospheric conditions. 

Energy Balance 


The 3 -step procedure for estimating T D , when Ed and/or skies were not 

cloudfree, is based on the energy balance equation, i.e. 

R n -G = H+LE (3) 

where R n is the net radiant flux density, G is soil heat flux density and H is sensible heat 
flux density (all in units of W m'~). The A£ (W m'~) is the latent heat flux density that is a 
product of the heat of vaporization X (J kg' ) and the rate of evaporation E (kg s' m" ). 

Eq. (3) neglects the horizontal advective flow of heat and water vapor and values of G, H 
and XE are positive when directed away from the surface. For a bare soil surface when 
E=0, H can be expressed as 


H = C v (t s -t a )/r a (4) 

I O 

where C v the volumetric heat capacity of air (J EC' m' ), t s is the soil surface temperature 
(EC), t a the air temperature (EC), and r a the aerodynamic resistance (s m' 1 ). 

Jackson et al. (1981) wrote the energy balance equation in terms of foliage-air 

temperature, 
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(tf-t a )=[r a (R n -G)/C v ] [(( 1 +r c /r a )/ {)+((! +r c /r a ) } ]- [ VPD/ {)+(( 1 +r c /r a ) } ] , 


( 5 ) 


where tf is the plant foliage temperature (EC), r c the canopy resistance (s m' 1 ) to vapor 
transport, ( the psychrometric constant (kPa EC' 1 ), ) the slope of the saturated vapor 
pressure-temperature relation (kPa EC' 1 ), and VPD the vapor pressure deficit of the air 
(kPa). For saturated bare soil, where r c =0 (the case of a free water surface), 

t S M!N= ta + [r a (R n -G)/ C v ] [(/ ()+()] -[VPD/ ()+()] , ( 6 ) 

and for dry bare soil, where r c =4 (analogous to complete stomatal closure), 

tsMAx - t a + [r a (R n -G)/C v ]. ( 7 ) 


where t S MiN is minimum soil surface temperature (EC) and t S MAx is maximum soil surface 
temperature (EC) for the given meteorological conditions at a given time of day. 

The on-site measurements necessary to solve Eqs. (6)-(7) are R n , VPD, t a and 
wind speed (Moran et al., 1994). If R n over bare soil is unavailable, it can be computed 
from other meteorological measurements (Jackson et al., 1981) and a value of G can be 
estimated as a function of R n , if necessary (Clothier et al., 1986). The value of r a is 
notoriously difficult to compute and can lead to high uncertainty in t S MiN and t S MAX 
(Stewart et al., 1994). However, for this application, since ET is being measured on site, 
it is possible to invert Eq. (7) for dates when ETd is near zero and empirically estimate r a 
for the site. This value of r a would theoretically be only dependent upon wind speed and 
would potentially be much more accurate than theoretical computation. 

With t sM iN and t S MAx computed from Eqs. (6) and (7) and measured soil surface 
temperature (Lmeas), Moran et al. (1994) determined that at midday 


E/Ep - (tsMAX -t s MEAs)/(tsMAX -t s MIN) 


( 8 ) 


2 

where E p is potential evaporation (W m’ ). Then, actual Ed can be estimated by 
multiplying the daily E p (E p d) by E/E p assuming E/E p = Ed/E p d , where 


Ed - (E/E p )E p d 


( 9 ) 


and 


E pD = (A Rn) + ((Cv VPD) /r a )/(A+y) (10) 

(Allen, 1986) with E p d converted to units of mm/day based on the latent heat of 

* s i 'j 'y i 

vaporization (2.45 10 J kg' ) and density of water (1.0 10' m kg' ). Finally, daily 
transpiration can be determined by 


Td— ETd - Ed. 


( 11 ) 


STUDY SITE, MATERIALS AND METHODS 
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The Soil Moisture Experiments 2004 (SMEX04) was conducted during the 
summer of 2004 in Arizona and Mexico to address overlapping science issues of the 
North American Monsoon Experiment (NAME) and soil moisture remote sensing 
programs. As part of SMEX04, two sites in the USDA ARS Walnut Gulch Experimental 
Watershed (WGEW) were instrumented with automated sensors to measure surface and 
atmospheric conditions. The WGEW is located in southeastern Arizona with a semiarid 
climate characterized by cool, dry winters and warm, wet summers. Mean annual 
precipitation is 356 mm and mean annual temperature is 17 EC. 

The Kendall site is dominated by herbaceous vegetation, predominately black 
grama ( Bouteloua eriopoda), sideoats grama ( Bouteloua curtipendula), three-awn 
(. Aristida sp.) and cane beardgrass ( Bothriochloa barbinodis). The soil at the site is 
dominately Stronghold (Coarse-loamy, mixed, thermic Ustollic Calciorthids) with slopes 
ranging from 4 to 9%. Only 9 km to the west, the Lucky Hills site is a shrub plant 
community dominated by creosotebush ( Larrea divaricata ), whitethorn Acacia ( Acacia 
constricta), mariola (Parthenium incanum), and tarbush ( Flourensia Cernua). The soil at 
this site is Luckyhills series (Coarse-loamy, mixed, thermic Ustochreptic Calciorthids) 
with 3 to 8% slopes (King et al., this issue; Skirvin et al., this issue). 

At each site, soil surface temperature was measured with an IRT at 5 -minute 
intervals and precipitation was measured at 1 -minute intervals. Volumetric soil moisture 
(0) was measured at 3 depths (5, 15 and 30 cm) with Vitel capacitance sensors at 5- 
minute intervals (Paige and Keefer, submitted; Keefer et al., this issue). Soil temperature 
was measured at 1-, 2-, 5-, 6-, 15- and 30-cm depths with thermocouples at 20-minute 
intervals. Meteorological data (including incoming solar radiation and soil heat flux) 
were measured at 5- and/or 20-minute intervals. These sites were also equipped with flux 
stations to measure evapotranspiration using a Bowen ratio technique at 20-minute 
intervals (Emmerich, 2003; Emmerich and Verdugo, this issue (b)). 

During the growing season in 2003, measurements of ET and T were made at the 
Lucky Hills shrub-dominated site. ET was monitored every twenty minutes using the 
Bowen ratio method (Emmerich, 2003), and shrub transpiration was measured every 
thirty minutes using the constant heat balance sapflow technique (Scott et al., in press). 
This shorter, but more comprehensive, data set was used to supplement and clarify the 
analysis of the 2004/2005 study at Kendall and Lucky Hills. 

Thus, data sets of ET, T (in the 2003 study), meteorological data, volumetric soil 
moisture at 5cm, surface temperature (from IRT), and soil temperatures at multiple 
depths were compiled to study the partitioning of E and T and the computation of plant 
WUE. These data were analyzed over an eighteen-month period in 2004 and 2005 
encompassing the dry/hot season, the North American monsoon and the dry/cool season, 
with particular attention to drying periods after storm events. Some analysis was 
conducted for the growing season only, defined as the time when plants were transpiring, 
which corresponds to a period from about DOY 220 to 280 at WGEW. Rainfall records 
for that period and the entire year in 2003, 2004 and 2005 are summarized in Table 1 . 

RESULTS 

As discussed in the previous section, Ia is theoretically related to both surface and 
atmospheric conditions. This sensitivity is illustrated by the response of Ia to a variety of 
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surface and atmospheric conditions associated with a spring storm at Kendall before the 
summer vegetation growth (Figure 1). A precipitation event on DOY 147 and 148 
resulted in a dramatic decrease in I a associated with an increase in ETd (due to increased 
soil moisture) and an associated decrease in available solar energy (due to cloudiness). 
For the clear-sky days that followed the storm event (DOY 149-152), Ia steadily 
increased as ETd decreased, finally reaching a value similar to that before the storm. 
However, cloudy conditions on the following day (DOY 153) resulted in another 
dramatic decrease in Ia without any significant change in soil moisture. This 
demonstrates the difficulty in interpretation of Ia, and introduces the rationale behind the 
approach used here. 

Apparent Thermal Inertia Related to Soil Moisture and Evanotransniration 

Results showed that Ia was not well related to soil moisture for representative 
summer, winter and spring storms in 2004 and 2005 at WGEW (Figure 2). At Kendall 
(grassland) and Lucky Hills (shrubland), the Ia decreased immediately with precipitation, 
but returned to its pre-storm value within days of the storm, depending on atmospheric 
conditions. In contrast, surface soil moisture (at 5 cm) reached a peak a day or two after 
the storm, but continued to decrease for weeks thereafter. 

For the same winter and spring storms (when transpiration was known to be zero 
because vegetation was senescent), the Ia related well with ETd (Figure 1 and Figure 3b). 
Generally, the Ia was inversely correlated with ETd and both measures returned to their 
pre-storm values within the same time period. For the 2004 summer storm (Figure 3a), 
the Ia post-storm recovery corresponded to a steep decline in ETd (related to Ed). This 
was followed by a more gradual decline in ETd, related to Td- This trend was confirmed 
by the Ia and Ed measurements made at Lucky Hills in 2003 (Figure 4). For two small 
summer storms, variation (decrease and recovery) in Ia corresponded directly to the 
measured increase and subsequent cessation of Ed- 

As one would expect from these results, the relation between soil moisture at 5 cm 
depth and ETd is weak in both winter and summer (Figure 5). The ETd is highly 
influenced by storm events and solar radiation, whereas the soil moisture has a less 
dramatic post-storm peak and steadily decreases until the next storm event. 

Partitioning E and T from ET with Apparent Thermal Inertia 

Based on the results in the previous sections, we postulated that the highest I A 
values were associated with cloudfree days when Ed was negligible. We also observed 
that Ia followed a seasonal trend in which higher values were obtained during the 
summer when solar radiation was at a maximum. To extract the days when Ed~0 and 
ETd~Td, it was first necessary to detrend the annual Ia time series. For ease of 
computation, this was done in three 6-month sets, as follows. 

For 6-month periods in 2004 and 2005, a polynomial was fit to the highest I A 
values in the data stream (Figure 6a). Then, an adjustment was made to all the values to 
remove the seasonal trend relative to the first Ia value in the data stream, resulting in a 
detrended Ia' (Figure 6b). Finally, a threshold was determined (10% of the highest value 
Ia' value) to select only the highest values of Ia' (Figure 6b). For dates which Ia' 
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exceeded that threshold, we presume that Ed was negligible and ETd~T d . Thus, To was 
estimated for selected dates for predominantly grassland (Kendall) and woody vegetation 
(Lucky Hills) over the time period 2004 and 2005 (Figures 6c and 6d). 

This process is premised on the assumption that the measurement of ts with the 
IRT at one site represents the conditions across the site, including both sunlit and shaded 
(under shrub cover) soils. Tuzet et al. (1997) described a “shade effect” that could affect 
the partitioning of ET due to variability of soil surface water availability. In 1990 at 
WGEW, electrical resistance sensors (ERS) soil moisture probes were deployed at Lucky 
Hills under three bare and three shrub-covered surfaces at 5-cm depths (Hymer et al., 
2000). The difference between surface soil moisture between bare and shrub-covered 
surfaces was only 0.2% (0.002 volumetric soil moisture) over the entire 1990 growing 
season (DOY 229-280). 

This process could be applied using near-surface soil temperature measurements, 
rather than IRT measurements of soil surface temperature. This might be preferable since 
the instrumentation is less expensive. We found that the amplitude of I A decreased with 
depth in the soil (Figure 7). Nonetheless, I A ' computed from soil temperatures at 1cm 
produced results similar to I A ' based on IRT measurements (compare results in Figure 8 
and Figure 6c). 

Deriving Tn when E ?© and/or Cloudy Sky Conditions 


To this point, T D has been derived from a combination of ET D measurements and 
diurnal surface temperature differences. This approach can be applied only when Ed~0 
and the sky is predominately clear over the daytime period. It was argued in a previous 
section that To could be determined for days when Ed?© and/or the sky conditions were 
not clear by using energy balance theory, with a site-calibrated aerodynamic resistance 
(r a ). The r a for Kendall and Lucky Hills was computed by inverting Eq. (7) to solve for r a 
using measurements of R n , G, t a and C v at 2:00 pm on clear sky days in June 2004 and 
2005 when ETd measurements were near zero. The results show that the r a value at both 
sites was relatively stable (ranging from 42 to 63 s m" 1 at Kendall and 41 to 68 s m’ 1 at 
Lucky Hills) over a range of wind speeds from 1 .7 to 6.6 m s' 1 (Figure 9). Similar results 
were found by Holifield Collins et al. (this issue) for Kendall using another dataset. The 
maximum r a values for all dates were used to solve the energy balance equations (6) and 
(7) to derive the maximum- and minimum-possible soil surface temperatures (Lmin and 
tsMAx) for midday (2:00 pm) at Kendall and Lucky Hills for input to Eq. (8) and 
determination of To (Eq. 11). 

An indirect test of the accuracy of the derived r a values is to compare the H 
derived by Eq. (7) with the H measured by the Bowen ratio instrumentation. For clear 
sky days during the growing season (DOY 220-280) at 2:00 pm in 2004 and 2005, the H 
derived from Eq. (7) based on empirical r a values and measured temperatures 
corresponded well with H measured independently by the Bowen ratio sensors with mean 
absolute differences (MAD) of 32 W m" 2 and 35 W m" 2 for Kendall and Lucky Hills, 
respectively (Figure 10). These MAD are comparable to differences found between 
values measured by two Bowen ratio systems at the same location (Houser et al., 1998). 

With this confidence in the theoretical approach, we computed t sM iN and t sMA x for 
each day during the growing seasons (DOY 220-280) at Kendall and Lucky Hills in 2004 
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and 2005 (e.g., Figure 1 1 a). Based on Eqs. (8)-(l 1), we estimated To for days when 
Ed ?© and/or skies were not cloudffee. We added a rule to this process to avoid the 
computationally possible, but illogical, situation in which Ed would decrease to zero 
during the few days after a storm and then increase to greater than zero without any 
further storm activity. The rule was that once Ed=0, it remained zero until the next 
rainfall event. Combined with the values of Td computed from (e.g. Figures 6c and 
6d), the result was a season-long estimate of Td and Ed with reasonable trends at both 
sites and both years (e.g. Figure 1 lb). 

Validation 2003 


Validation of this approach was based on measurements from a study by Scott et 
al. (in press) designed to measure ecosystem Td as well as EToat Lucky Hills in 2003. 
The same series of steps illustrated in Figure 6 for Kendall 2004 were followed to 
estimate T for clear-sky days when E=0 using Ia’ at Lucky Hills in 2003. For all other 
days, T was computed using the energy balance equations (4)-(l 1) as described in 
Section x. The correlation between T estimated with Ia’ and T measured with a sapflow 
technique was good (illustrated with asterisks in Figure 12, MAD = 0.29 mm/d). The 
correlation between T estimated with the semi-empirical energy balance approach was 
not as successful (illustrated with circles in Figure 12, MAD = 0.64mm/d). The worst 
result (absolute difference of 2.1 mm/d) was obtained when the day was cloudy for only 
part of the day (DOY 236 in Figure 12b). Under these conditions, the surface 
temperatures computed at 5 am and 2 pm (Lmin and Lmax) did not necessarily 
characterize the Ed and the computation of E p d was unreliable. 

The ratio of T and ET for the growing season (Ts/ETs) in 2003 (DOY 220 to 280) 
based on measurements by Scott et al. (in press) was 0.84 (Figure 13). The Ts/ETs based 
on the estimation approach presented here was 0.76, which was within 10% of the 
measured value. In both cases, T s was the dominant component of ET S during the North 
American monsoon season at WGEW. Preliminary estimates of Ts/ETs for the growing 
seasons in 2004 and 2005 show that Ts/ETs ranged from 0.50 to 0.79 (Figure 14) and the 
value depended on rainfall and meteorological conditions that characterized the season 
(Table 1). Based on analysis of 2004 and 2005 results, Ed is close to zero within a few 
days of the storm event during the growing season. 

CONCLUSION 

The main contribution of this work was a new method to partition on-site 
measurements of ETd into daily E and T based on the low-cost addition of an IRT to 
existing eddy covariance and Bowen ratio stations. We showed that reasonable estimates 
of To were obtained for clear-sky days when Ed was found to be zero based on the 
magnitude of the apparent thermal inertia (Ia). This finding is the most valuable and 
reliable aspect of this approach. The approach is based on instrumentation that can be 
maintained in place continuously for years with no more expertise and effort than is 
already required for deployment of energy flux stations. 

A theoretical supplement was added to compute Td on days when the U-based 
approach is not applicable (e.g. cloudy days or when Ed t©), resulting in season-long 


9 



estimates of T D and Ed- The theoretical energy balance approach should be refined to 
account for days with variable cloud conditions. Alternatively, it may be possible to use 
some empirical information gained from days when Ed is known to be zero to derive an 
empirical relation between To and another site measurement (e.g. surface (5 cm) soil 
moisture) to partition ETd on those days. 

Future work should continue to explore the soil water availability in grass- and 
shrub-dominated ecosystems like Kendall and Lucky Hills to better understand the 
potential impact of woody plant encroachment in grasslands across the Western U.S. 
There is a great deal of evidence that the Ts/ETs ratio is sensitive not only to woody plant 
cover, but to the variation in amounts, frequency and timing of rainfall (Reynolds et al., 
2000; Guswa et al., 2002; Bhark and Small, 2003; Loik et al., 2004). These preliminary 
results will be combined with results from similar data collected in 2006 (a relatively wet 
monsoon season) to further study this link between Td/ETd and such variations in rainfall 
at WGEW. 
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Figure and Table Captions: 

Table 1 . Rainfall summary (mm) for the growing season at WGEW during years 2003- 
2005 at Kendall and Lucky Hills sites. 

Figure 1. Comparison of volumetric soil moisture at 5 cm at 2:00 pm (0, %), the 
apparent thermal inertia derived from IRT measurements at 2:00 pm and 5:00 am (Ia, °C, 
Eq. 2), and daily ET measured with the Bowen ratio method (mm/d, multiplied by 10 for 
presentation) for a storm event at Kendall in 2005. Circles indicate the cloudffee days and 
bars represent daily precipitation (mm). 

Figure 2. Comparison of Ia (Eq. 2) with volumetric soil moisture (0) at 5cm at 2:00 pm 
for summer and winter storms at Kendall, followed by a long series of cloudffee days. 
Similar results were found (though not shown here) for Lucky Hills and for the spring 
storm in 2005. Bars represent daily precipitation (mm). 

Figure 3. Comparison of Ia (Eq. 2) with daily ET (multiplied by 10 for presentation) for 
summer and winter at Kendall (the spring storm was presented in Figure 1), followed by 
a series of cloudffee days. Similar results were found (though not shown here) for Lucky 
Hills. Bars represent daily precipitation (mm). 

Figure 4. Comparison of Ia (Eq. 2) with daily E (multiplied by 10 for presentation) for 
two storms in 2003 at Kendall, when Ed was determined ffom the difference between 
ETd (using Bowen ratio) and To (using sapflow technique). Bars represent daily 

precipitation (mm). 

Figure 5. The weak relation between surface soil moisture (at 5 cm) and daily ET D 
(multiplied by 10 for presentation) for a series of storms in the dry/hot season (DOY120- 
180), the North American monsoon (DOY 200-270), the growing season (DOY 220-280) 
and the dry/cool season (near DOY 360). Bars represent daily precipitation (mm). 

Figure 6. An illustration of the steps taken to partition ET using Ia at Kendall in 2004. 

(a) A polynomial was fit to the highest Ia vales, (b) A threshold was set to discriminate 
the highest detrended Ia values (Ia')- (c) For dates when Ia' exceeded the threshold, then 
ETd=Td. (d) Values of To for Lucky Hills were derived using the same process 
illustrated at Kendall in Figures 6a-6c. Bars represent daily precipitation (mm). 

Figure 7. I A computed ffom soil temperature measurements at the surface (solid line), 
and at depths of 1 cm, 2 cm and 6 cm in the soil (with thermocouples). Similar results 
were found (though not shown here) for Lucky Hills and other storms. The bars represent 
daily precipitation (mm). 

Figure 8. Daily transpiration at Kendall in 2004 derived ffom Ia' using soil temperature 
at 1 cm instead of IRT measurements, following the steps illustrated in Figure 6a-6c. 
These results can be compared to results presented in Figure 6c using the IRT. 
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Figure 9. The impact of wind speed (U) on estimates of aerodynamic resistance (r a ) at 
Kendall and Lucky Hills , where r a was computed by inverting Eq. (7) with measurements 
made at 2:00 pm on clear sky days in June 2004 and 2005when ETd was near zero. 

Figure 10. Comparison of sensible heat flux (H) at 2:00 pm at Kendall and Lucky Hills 
measured by the Bowen ratio instrumentation on clear-sky days in 2004 and 2005 during 
the growing season (DOY 220 to 280) and “estimated” H derived from Eq. (4) based on 
on-site measurements and the empirically derived r a (from Figure 9). The MAD is the 
mean absolute difference between measured and estimated values. 

Figure 11. Values of a) t s , t S MiN and Lmax and b) T and E for each day during the 
growing season (DOY 220-280) at Lucky Hills in 2005. Similar results were found, 
though not presented here, for Kendall 2004 and 2005 and Lucky Hills 2004. The bars 
represent daily precipitation (mm). 

Figure 12. (a) Comparison of transpiration (T) estimated with this approach and T 
measured with a conventional sapflow technique (from Scott et al., in press) and (b) the 
difference between measured and estimated T (T meas -Testim)- The results estimated with 
Ia’ are illustrated with an asterisk, and results estimated with the semi-empirical energy 
balance approach are illustrated with a solid circle. The bars represented daily 
precipitation (mm). 

Figure 13. Ratio of T and ET for the growing season (DOY 220 to 280) at Lucky Hills 
in 2003, based on measurements (Meas) and on the proposed estimation approach 
(Estim). 

Figure 14. Preliminary estimates of the ratio of season-long transpiration and 
evapotranspiration (Ts/ETs) for the growing season (DOY 220 to 280) at Kendall and 
Lucky Hills sites in (a) 2004 and (b) 2005. 
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Table T Rainfall summary (mm) for the growing season at 
WGEW during years 2003-2005 at Kendall and Lucky Hills 

sites. 

Site and Year 

Annual 

Growing Season DOY 
220-280 

Total 

Precipitation 

Total 

precipitation 

# of storms 

Lucky Hills 2003 

245.9 

93.2 

15 

Kendall 2004 

294.9 

67.7 

12 

Lucky Hills 2004 

219.4 

58.4 

12 

Kendall 2005 

162.3 

61.0 

16 

Lucky Hills 2005 

223.2 

87.2 

13 




Kendall 2005 Spring Storm 



DOY 


Figure 1. Comparison of volumetric soil moisture at 5 cm at 2:00 pm (0, %), the apparent thermal inertia 
derived from IRT measurements at 2:00 pm and 5:00 am (I A , °C, Eq. 2), and daily ET measured with the 
Bowen ratio method (mm/d, multiplied by 10 for presentation) for a storm event at Kendall in 2005. Circles 
indicate the cloudfree days and bars represent daily precipitation (mm). 



Figure 2. Comparison of I A (Eq. 2) with volumetric soil moisture (0) at 5cm at 2:00 pm for summer and 
winter storms at Kendall, followed by a long series of cloudfree days. Similar results were found (though 
not shown here) for Lucky Hills and for the spring storm in 2005. Bars represent daily precipitation (mm). 
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Kendall 2004 Winter Storm 
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Figure 3. Comparison of I A (Eq. 2) with daily ET for summer and winter at Kendall (the spring storm was 
presented in Figure 1), followed by a series of cloudffee days. Similar results were found (though not 
shown here) for Lucky Hills. Bars represent daily precipitation (mm). 
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Lucky Hills 2003 
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Figure 4. Comparison of I A (Eq. 2) with daily E for two storms in 2003 at Kendall, when E D was 
determined from the difference between ET D (using Bowen ratio) and T D (using sapflow technique). Bars 
represent daily precipitation (mm). 
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Figure 5. The weak relation between surface soil moisture (at 5 cm) and daily ET D (multiplied by 10 for 
presentation) for a series of storms in the dry/hot season (DOY 120- 180), the North American monsoon 
(DOY 200-270), the growing season (DOY 220-280) and the dry/cool season (near DOY 360). Bars 
represent daily precipitation (mm). 
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Figure 6. An illustration of the steps taken to partition ET using I A at Kendall in 2004. (a) A polynomial 
was fit to the highest I A vales, (b) A threshold was set to discriminate the highest detrended I A values (I A ). 
(c) For dates when I A exceeded the threshold, then ET D =T D . (d) Values of T D for Lucky Hills were derived 
using the same process illustrated at Kendall in Figures 6a-6c. Bars represent daily precipitation (mm). 
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Kendall 2004 Summer Storm 



Figure 7. I A computed from soil temperature measurements at the surface (solid line), and at depths of 1 
cm, 2 cm and 6 cm in the soil (with thermocouples). Similar results were found (though not shown here) for 
Lucky Hills and other storms. The bars represent daily precipitation (mm). 


Kendall 2004 (using t s at 1 cm) 



DOY 


Figure 8. Daily transpiration at Kendall in 2004 derived from I A ' using soil temperature at 1 cm instead of 
IRT measurements, following the steps illustrated in Figure 6a-6c. These results can be compared to results 
presented in Figure 6c using the IRT. 
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Kendall 2004 2005 June 



Lucky Hills 2004 2005 June 



Figure 9. The impact of wind speed (U) on estimates of aerodynamic resistance (r a ) at Kendall and Lucky 
Hills , where r a was computed by inverting Eq. (7) with measurements made at 2:00 pm on clear sky days in 
June 2004 and 2005when ET D was near zero. 
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Kendall 2004 2005 



Lucky Hills 2004 2005 



Figure 10. Comparison of sensible heat flux (H) at 2:00 pm at Kendall and Lucky Hills measured by the 
Bowen ratio instrumentation on clear-sky days in 2004 and 2005 during the growing season (DOY 220 to 
280) and “estimated” H derived from Eq. (4) based on on-site measurements and the empirically derived r a 
(from Figure 9). The MAD is the mean absolute difference between measured and estimated values. 
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Lucky Hills 2005 Growing Season 
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Figure 11. Values of a) t s , t sM iN and t sM Ax and b) T and E for each day during the growing season (DOY 
220-280) at Lucky Hills in 2005. Similar results were found, though not presented here, for Kendall 2004 
and 2005 and Lucky Hills 2004. The bars represent daily precipitation (mm). 
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Lucky Hills 2003 Growing Season 
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Figure 12. (a) Comparison of transpiration (T) estimated with this approach and T measured with a 
conventional sapflow technique (from Scott et al., in press) and (b) the difference between measured and 
estimated T (T meas -T estim ). The results estimated with I A ’ are illustrated with an asterisk, and results 
estimated with the semi-empirical energy balance approach are illustrated with a solid circle. The bars 
represented daily precipitation (mm). 
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Figure 13. Ratio of T and ET for the growing season (DOY 220 to 280) at Lucky Hills in 2003, based on 
measurements (Meas) and on the proposed estimation approach (Estim). 



Figure 14. Preliminary estimates of the ratio of season-long transpiration and evapotranspiration (T S /ET S ) 
for the growing season (DOY 220 to 280) at Kendall and Lucky Hills sites in (a) 2004 and (b) 2005. 
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